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We investigate the topology of the spin-polarized charge density in bcc and fee iron. While the 
total spin-density is found to possess the topology of the non-magnetic prototypical structures, in 
some cases the spin-polarized densities are characterized by unique topologies; for example, the 
spin-polarized charge densities of bcc and high-spin fee iron are atypical of any known for non- 
magnetic materials. In these cases, the two spin-densities are correlated: the spin-minority electrons 
have directional bond paths with deep minima in the minority density, while the spin-majority 
electrons fill these holes, reducing bond directionality. The presence of two distinct spin topologies 
suggests that a well-known magnetic phase transition in iron can be fruitfully reexamined in light 
of these topological changes. We show that the two phase changes seen in fee iron (paramagnetic 
to low-spin and low-spin to high-spin) are different. The former follows the Landau symmetry- 
breaking paradigm and proceeds without a topological transformation, while the latter also involves 
a topological catastrophe. 

PACS numbers: 75.50.Bb, 75.30.Kz 



Bader's topological theory of molecular structure, 
Atoms in Molecules (AIM) , has been successfully applied 
to a variety of crystalline systems [IJIllllllllElHIHJISI 
[inilll]. It has been employed to investigate the nature of 
bonding in materials ranging from high temperature al- 
loys to biological systems. These have yielded surprising 
results, such as second neighbor bond paths in B2 ionic 
crystals [5] and transition metal aluminides U, with the 
magnitude of the latter correlating to failure properties 
Other studies have used bond path properties to of- 
fer first principles explanations of stress-induced failure 
in brittle and ductile alloys [7], as well as shear elastic 
constants in a variety of pure metals and alloys [H [9l [10] . 
Building on these ideas and using the rigorous definitions 
of bond paths afforded by the theory, the anomalous be- 
havior of iridium under shear was also explained . 

One of the attractive features of AIM is its reliance 
on the charge density, a quantum mechanical observable, 
that is most often calculated but can, in principle, be 
measured via X-ray diffraction techniques [12l[T3 l [T4 l ll5j. 
In a similar fashion, the spin-polarized charge density is 
an observable that can be calculated or measured using 
spin-polarized neutron diffraction. Despite the informa- 
tion and insights that have come from topological investi- 
gations of the total charge density, the same analysis has 
yet to be performed on spin-polarized densities. Here, we 
report the results from the first such studies, exploring 
the spin-minority and spin-majority topologies of body- 
centered-cubic (bcc) and face-centered-cubic (fee) iron. 

This first application of AIM to spin-density sheds light 
on the origins of the magnetic phase transitions of fee 
iron. It is argued that this system undergoes two distinct 
phase transitions during volume expansion: a second- 
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order phase change, from a paramagnetic to a low-spin 
state, occurs at smaller volumes. It is coincident with a 
change in the charge density at the critical points without 
a topological transformation. At larger volumes, the low- 
spin state changes to high-spin through a discontinuity 
in the moment, where the topology also transforms. 

The framework of AIM can be applied in the same 
manner to both the total and spin-polarized densities. 
It is known from the Hohenberg-Kohn theorem that all 
ground-state molecular properties are a consequence of 
its charge density p{f) [TB], a scalar field in three spa- 
tial dimensions. Bader noted p{'r) must also contain the 
essence of a molecule's structure, which can be described 
topologically. 

The topology of a general scalar field can be charac- 
terized in terms of its critical points (CPs) , the zeros of 
the gradient of this field. There are four kinds of CP 
in a three-dimensional space: a local minimum, a local 
maximum and two kinds of saddle points. These CPs are 
conventionally denoted by an index, which is the num- 
ber of positive curvatures minus the number of negative 
curvatures. For example, a minimum CP has positive 
curvature in three orthogonal directions; therefore it is 
called a (3, 3) CP. The first number is simply the num- 
ber of dimensions of the space, and the second number 
is the net number of positive curvatures. A maximum is 
denoted by (3, -3), since all three curvatures are negative. 
A saddle point with two of the three curvatures negative 
is denoted (3, -1), while the other saddle point is a (3, 1) 
CP. 

Through extensive studies of molecules and crystals, 
Bader and Zou and Bader showed that it was possible to 
correlate topological properties of the charge density with 
elements of molecular structure and bonding [U |2]. In 
particular, a bond path was shown to correlate with the 
ridge of maximum charge density connecting two nuclei. 
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such that the density along this path is a maximum with 
respect to any neighboring path. The existence of such 
a ridge is guaranteed by the presence of a (3, -1) CP 
between bound nuclei. As such, this CP is sometimes 
referred to as a bond CP. 

Other types of CPs have been correlated with other 
features of molecular structure. A (3, 1) CP is topo- 
logically required at the center of ring structures, e.g. 
benzene. Accordingly, it is designated a ring CP. Cage 
structures are characterized by a single (3, 3) CP and 
again are given the descriptive name of cage CPs. The 
locations of the atomic nuclei always coincide with a max- 
imum, a (3, -3) CP. Hence it is conventionally called an 
atom CP. 

In this study, iron's spin-polarized charge densities 
were calculated using the Vienna ab-initio simulation 
package (VASP) version 4.6 [171 HI]) which is based on 
the projector augmented wave method (PAW) |T9l l20j . 
The Perdew-Wang (PW91) generalized gradient cor- 
rections [51] and the Vosko-Wilk-Nusair interpolation 
[22j [23| of the correlation part of the exchange-correlation 
functional were included. 

Calculations performed on non-magnetic, monotonic, 
bcc transition metals show that the ground state total 
charge densities share the same topology as the prototype 
structure, W. Each atom has eight bonds to the nearest 
neighbor atoms, with no known case of second neighbor 
bond paths in elemental bcc materials. This also holds 
true for the total charge density of ferromagnetic bcc (a) 
Fe; however it does not hold for the spin-density topology. 

Two phases of bcc Fe with lattice constants (moments) 
of 2.87 A(2.22 hb) H and 2.83 1(2.21 /xb), re- 
spectively, were investigated, yielding identical topolo- 
gies. The total charge density topology is prototypical 
bcc, as can be seen in Fig. [iJV, where a contour plot 
through second neighbor atoms in a (100) plane is shown. 
A cage point, circled in Fig.[lJV, is present. Similarly, the 
spin-minority density reflects the topology of the total. 
Fig. [l]3. When the spin-majority density is examined, 
however, a different picture emerges. Now, instead of a 
cage point between second neighbor atoms, a bond point 
is found, Fig. [Tp, giving rise to 14 bond paths in the 
spin-majority topology. 

Our results are consistent with the band theory of fer- 
romagnetism, wherein, the d-band splits into bonding, 
spin-minority, and less bonding, spin-majority, bands. As 
bcc Fe has Oh symmetry about the atoms, the d-orbitals 
reduce as T2g and Eg at the band F-point. The spin- 
minority electrons predominately populate the triply de- 
generate state, to produce 8 (T-bonds directed towards 
cube corners. The electrons in the spin-majority band 
are more equitably distributed between the T2g and Eg 
bands. The inclusion of electrons of Eg symmetry in- 
creases the charge density at the cube faces, e.g., second 
neighbor bond paths are formed. These are, however, 
weakly directional. 

Directionality has been quantified by two components, 
using the quadratic surface constructed from the Hessian 



of the charge density (or spin-polarized charge density) at 
the bond CP [H [9l [101 IH] ■ An elliptic cone surrounding 
the bond path is a quadratic surface consistent with the 
signs of the principal curvatures at the bond CP. The ex- 
treme angles of the cone with respect to the plane normal 
to its axis are given by 

tana = , r^-^-^ and tan/3 — . r^-^ (i) 

V p\\\\ V p\\\\ 

where p|||| is the principal curvature parallel to the bond 
path and p±± and p±'±> are the two principal curvatures 
perpendicular to the bond path, which are degenerate in 
the bcc structure. This notation is employed to refiect 
the fact that the terms are all second derivatives of p{r) in 
the direction indicated by the subscripts. The definition 
in Eq. [l] is trivially generalized to directionality for the 
spin-polarized charge density by affixing a spin index. 

This definition was motivated by the desire to associate 
a "distance" to bond breaking, i.e. the vanishing of a 
bond CP and the transformation of the elliptic cone to 
intersecting planes. This will occur when one (or both) of 
the principal curvatures perpendicular to the bond path 
vanishes. By definition, then, bonds with larger values of 
directionality are further from a topological instability. 
The degenerate components of the spin-minority bond 
directionality are found to be 0.565, while those of the 
spin-majority are 0.472 and 0.227 for the first and second 
neighbor bond paths, respectively. In comparison, the 
bond directionality is 0.597 in W and 0.502 in bcc Fe. 

When non-magnetic fee metals were examined, only 
first neighbor bond were found in the total charge den- 
sity, Fig.[2j Each atom has a total of 12 bond paths, and 
both the octahedral and tetrahedral holes contain a cage 
point. However, as with magnetic bcc Fe, the spin densi- 
ties give rise to distinct topologies. In the case of fee (7) 
Fe, three magnetic phases are observed, high-spin, low- 
spin, and paramagnetic (HS, LS, and PM respectively) 
phases, each characterized by a different range of lattice 
constants [MUSZl UHl 1301131] • From the VASP calcu- 
lations we find the HS phase corresponds to a > 3.564 A, 
while the LS, low volume phase, occurs for a < 3.563 A, 
with the PM phases equilibrating at a = 3.45 A. 

The LS/HS phase change proceeds through a magne- 
tovolume instability [28l [29]. As the lattice parameter 
is increased from its LS phase (a = 3.465 A) value, the 
moment increases smoothly from 0.75 /is to 1.3 fj,B at 
a — 3.563 A. Then over a 0.001 A increase in lattice 
constant a near doubling of the moment to 2.5 fiB at 
a = 3.564 A results. Our results agree well with other 
computational [301 131] and experimental studies [321 133] ■ 
Less clear is the PM/LS phase change. The situation 
yields ambiguous results if p is used as an order param- 
eter, with the transition appearing as first [SS] [55] or 
second [30] order. Small changes in total energy associ- 
ated with this phase change, 0.001 eV/atom pOl, make 
the correct assignment of moment problematic. 

The topology of the PM phase is identical to the fee 
prototype, where a and (3 are not equivalent. In the fee 




FIG. 1: Total (A), spin-minority (B), and spin-majority (C) charge density contour plots in bcc Fe. The cut-plane passes 
through second neighbor atoms. Darker regions correspond to low charge density. Cage points, circled in white, are present in 
the total (A) and spin minority (B) charge densities. Second neighbor bond points (one circled in white) replace these in the 
spin-majority (C) charge density. 




FIG. 2: Topology of fee prototype. Only bond points (white) 
and maxima (gray) are shown for clarity. Paramagnetic, LS 
and the total and spin-majority density in HS fee Fe display 
this topology. (Though not shown, cage points are present in 
the tetrahedral holes.) 



structure, the directions of principal curvature perpen- 
dicular to the bond path at the bond CP point toward 
the octahedral and tetrahedral holes are denoted as poo 
and ptt respectively. Now, the components of direction- 
ality become 



tan a = 



and tan /3 = 




(2) 



For Fe and Cu respectively, the values of tan a (tan (3) 
are found to be 0.562 (0.456) and 0.463 (0.396). This in- 
crease is expected, as directionality decreases in the late 
transition metals [HI [9l [TOl [11] , which corresponds to the 
filling of anti-bonding orbitals. Similarly the topology of 
the LS phase is identical to the prototype. While the 
transition from PM to LS shows no topological change, 
the charge density at the bond CPs changes rapidly. Both 
the PM and LS phases of fee Fe have bonds with direc- 
tionality greater than that of bcc Fe and prototypical fee 



topologies. After a transformation to the HS phase, a 
different picture emerges. 

The topology of the spin-majority and total charge 
density in the HS phase is identical to the PM and LS 
phases. Bond directionality is, however, reduced, in ac- 
cordance with the band theory of ferromagnetism. Both 
the total and spin-majority have bond paths with direc- 
tionality nearly equal to that of bcc Fe, with tan a (tan /3) 
= 0.522 (0.337) and tana^ (tan/jT) = 0.466 (0.402), re- 
spectively. This strong similarity is a result of the filling 
of both T2g and Eg, orbitals by the spin-majority, as in 
the bcc metal. 

Interestingly the population of the d-orbitals reducing 
as Eg in the spin-majority band results in a change in 
the spin-minority topology. Here, non-nuclear maxima 
form in the tetrahedral holes. In Fig. [3] these can be seen 
acting to form a pseudo-bcc topology by allowing each Fe 
atom to form 8 cr-bonds. This is due to the fact that d- 
orbitals reducing as T2g vastly outnumber those reducing 
as Eg at moments in excess of 2.5 ps- As such, 12 bonds 
cannot be formed, as that would require both the T2g and 
Eg representations. In order to stabilize the structure, 8 
directional bonds, tana = 1.124, are created. 

This change in topology is responsible for the magne- 
tovolume instability. Only two stable topologies exist in 
fee Fe: one corresponds to that observed in the PM-LS 
phase, while the other is that of the spin-minority band 
of the HS phase. The PM/LS (prototypical fee) topol- 
ogy can be formed so long as /i < l.'ip.B- Moments in 
this regime still allow the spin-minority to use the Eg 
representation for directional bond path formation. In 
order for this topology to transform to the HS topology, 
a CP annihilation-a topological catastrophe-must occur 
within the spin- minority topology. The spin-majority 
topology, on the other hand, remains unchanged through 
the phase transition. 

The transformation of the spin-minority topology re- 
quires the 12 bond points about each atom in the LS 
phase to become ring points. In conjunction, the tetrahe- 
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FIG. 3: Topology of the minority spin-density in HS fee Fe. 
Only bond points (white) and maxima (gray) are shown for 
clarity. The appearance of non-nuclear maxima allows the Fe 
atoms to adopt a bcc-like topology. These replace the cage 
points found in the tetrahedral holes of the prototype, and 
the bond points replace the ring points of the prototype. The 
inset shows the tetrahedron of bond points around a non- 
nuclear maximum. 

dral cage points must transform to non-nuclear maxima. 
At the same time, ring points become the bond points 
shown in Fig. [3] During the catastrophe, the curvature 
at the bond points will become flat in one direction, e.g. 
a zero eigenvalue of the Hessian (and a vanishing of the 
Gaussian curvature). Thus, when the population in the 
Eg representation drops (1.3/iB < /i < 2.5/iB), bond 
paths in the spin-minority band are weakened and the 
resultant topology is unstable. 

By way of example, the directionality of the bond 
points in the minority spin-density, measured from the 
tetrahedral hole, drops from tan 13^ = 0.380 to 0.240, 
as the moment is changed from 1.3 hb to 2.1 jj.^ at 



a = 3.563A. If, however, the moment is allowed to in- 
crease to 2.5 /iB, the number of bond paths is reduced 
to 8, but their directionality is increased to 1.20. Thus, 
any moment in the range of l.SfJ-B < /i < 2.5/^b will not 
be observed. Instead, the topology will spontaneously 
relax, resulting in a discontinuous change in moment. In 
a one-electron picture, it is the shift of electrons from the 
spin-minority to the spin-majority Eg bands that brings 
about the phase transformation. 

In summary, we have applied Bader's AIM theory to 
the spin-density in bcc and fee iron. It has been shown 
that the spin-majority and spin-minority topologies are 
different in both bcc Fe and the HS phase of fee Fe. In 
each case, the two topologies accumulate charge in differ- 
ent regions due to spin-spin correlation. While the spin- 
minority density in bcc Fe has cage points between sec- 
ond neighbor atoms, the spin-majority density has bond 
points. Here it is the spin-minority topology that resem- 
bles the total. In HS fee Fe this picture is reversed. The 
spin-majority topology is typical of fee metals with bonds 
points lying between the 12 nearest neighbor atoms and 
cage points in both the tetrahedral and octahedral holes. 
The spin- minority topology, however, has cage points 
in only the octahedral holes. Non-nuclear maxima are 
present in the tetrahedral holes. These serve to form a 
pseudo-bcc topology with 8 bond paths per atom. 

Using the fact that the topology can only transform 
through a catastrophe allows us to discern the two mag- 
netic phase changes in fee Fe. While describing these 
using fj, as a, local order parameter has offered conflicting 
results computationally, the topological treatment is un- 
ambiguous. A topological phase change does not occur 
in the PM/LS phase change. The LS/HS transition is, 
however, accompanied by a catastrophe, and the phase 
change can be conveniently followed by examining the 
abrupt change in the Bader bonding graph topology. We 
are grateful to the Office of Naval Research for financial 
support of this work. 
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